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The B°-B° oscillation frequency has been measured with a sample of 23 million BB pairs collected 
with the BABAR detector at the PEP-II asymmetric B Factory at SLAC. I 11 this sample, we select 
events in which both B mesons decay semileptonically and use the charge of the leptons to identify 
the flavor of each B meson. A simultaneous fit to the decay time difference distributions for opposite- 
and same-sign dilepton events gives Amj = 0.493 ± 0.012 (stat) ± 0.009 (syst) ps _1 . 

PACS numbers: 13.25.Hw, 12.39.Hg 


A precise measurement of the B°-B° oscillation fre¬ 
quency A md, is of fundamental importance as a direct 
measure of the Cabibbo-Kobayashi-Maskawa (CKM) ma¬ 
trix element \Vtd\ IQ- When combined with measure¬ 


ments of the B s -B s oscillation frequency, it provides 
a stringent constraint on the Unitarity Triangle of the 
CKM matrix. 

In this Letter, we present a measurement of the time 


4 


dependence of B°-B a mixing using data collected with 
the BABAR detector at the PEP-II asymmetric energy 
e + e _ collider operated at or near the T(4S) resonance. 
The data sample, recorded in the years 1999-2000, cor¬ 
responds to an integrated luminosity of 20.7fb~ 1 on the 
T(4S') resonance ( on-resonance sample), and 2.6 fb -1 col¬ 
lected 40 MeV below the T(4S) resonance ( off-resonance 
sample). BB pairs from the T(4S) decay move along the 
high-energy beam direction (z) with a nominal Lorentz 
boost (/Ty) = 0.55. Therefore, the two B decay vertices 
are separated by about 260 /im on average. 

The measurement technique is based on the identifica¬ 
tion of events containing two leptons from semileptonic 
decays of B mesons. The flavor of the B mesons at the 
time of their decay is determined or “tagged” by the 
charge of the leptons. Thus, for T(4S) resonance decays 
into B°B° pairs, neglecting backgrounds, opposite-sign 
(+) and same-sign (—) lepton pairs correspond to un¬ 
mixed and mixed events, respectively. Because the B°B° 
pair is in a coherent P- wave state, the time evolution of 
the B mesons is a function of the proper time difference 
At between the two B decays: 

e -|At|/r 

S±(At\Amd) = -(1 ± cos ArridAt), 

4 T 

where r is the B° lifetime (AT = 0 is assumed). The 
corresponding time-dependent asymmetry is (5+(At) — 
<S_(Af))/(<S+(At) + S-(At)) = cos ArridAt. 

This simple picture is modified by the effects of de¬ 
tector resolution and the presence of backgrounds. The 
most important background, about 50% of the sample, 
is due to B + B~ events, which are not removed by the 
event selection critieria. The fraction of B + B~ events is 
determined from the data itself in order to reduce sys¬ 
tematic uncertainties. Other non-negligible backgrounds 
are leptons from the b —> c —» t decay chain ( cascade 
decays), which are also the main source of wrong tags, 
and hadrons that are misidentified as leptons. Signal 
and background probability density functions (PDF) for 
opposite- and same-sign events are included as additional 
terms in the full PDF. The corresponding likelihood func¬ 
tion, combining opposite- and same-sign dilepton events, 
is maximized to determine Anq. 

The BABAR detector is described in detail elsewhere [|j . 
Charged particle tracking is provided by a 5-layer, 
double-sided silicon vertex tracker (SVT) and a 40- 
layer drift chamber (DCH), both operating inside a 1.5- 
T super-conducting solenoidal magnet. The CsI(Tl) 
electromagnetic calorimeter (EMC) detects photons and 
electrons. Particle identification is provided by a ring¬ 
imaging Cherenkov detector (DIRC) and specific ioniza¬ 
tion measurements dE/dx in the DCH. Muons are identi¬ 
fied with the instrumented flux return (IFR), segmented 
to contain resistive plate chambers. 

Events are selected by requiring more than 5 recon¬ 
structed charged tracks, at least 3 of which must originate 


from the interaction region and be reconstructed in the 
DCH. Continuum the normalized second Fox-Wolfram 
moment p] be less than 0.4 and the event aplanarity be 
greater than 0.01. Two-photon events, as well as resid¬ 
ual radiative Bhabha events with a large amount of miss¬ 
ing energy, are rejected by requiring the invariant mass 
squared of the event to be greater than 20GeV 2 /c 4 . 

Electrons are selected by requirements on the ratio of 
the energy deposited in the EMC to the momentum mea¬ 
sured in the DCH, the lateral shape of the energy depo¬ 
sition in the EMC, and dE/dx in the DCH. Muons are 
identified on the basis of the energy in the EMC, as well 
as the strip multiplicity, track continuity and penetration 
depth in the IFR. Lepton candidates consistent with the 
kaon hypothesis as measured in the DIRC are rejected. 
Electron (muon) selection efficiencies and misidentifica- 
tion rates at high momentum are about 92% (75%) and 
0.15% (3%), respectively. 

Photon conversions are rejected by combining each 
electron candidate with all other oppositely-charged elec¬ 
tron candidates in the event, selected with looser crite¬ 
ria, and applying requirements on the invariant mass and 
distance of closest approach of the pair in the transverse 
plane and along the beam direction. Leptons from J/if> 
and ij)(2S) decays are identified by pairing them with all 
other oppositely-charged candidates of the same lepton 
species, selected with looser criteria, and rejecting the 
whole event if any combination has an invariant mass 
within the J/ip or ip(2S) mass regions. 

Events with at least two leptons are retained and the 
two highest momentum leptons in the T(45 l ) rest frame 
are used in the following. 

The two lepton tracks and a beam spot constraint are 
used in a vertex fit to find the primary vertex of the event 
in the transverse plane. The positions of closest approach 
of the two tracks to this vertex in the transverse plane 
are computed and their z coordinates are denoted as z\ 
and Z‘ 2 , where the subscripts 1 and 2 refer to the highest 
and second highest momentum leptons in the T(4S) rest 
frame. The difference A z is defined as Az = 21 — 22 - The 
time difference At is computed from Az and the nominal 
boost as At = Az/{/3"/)c. 

The modeling of the resolution function 1Z is a cru¬ 
cial element of the A irid measurement. To improve the 
Az (and therefore At) resolution, reduce the fraction of 
incorrectly measured tracks, and minimize related sys¬ 
tematic uncertainties, charged tracks are required to sat¬ 
isfy the following criteria. Lepton candidates must have 
a distance of closest approach (DOCA) with respect to 
the nominal beam position of less than 1 cm in the trans¬ 
verse plane and less than 6 cm along the beam direction, 
at least 12 hits in the DCH, at least four 2 -coordinate 
hits in the SVT, a momentum in the range between 0.7 
and 2.5GeV/c in the T(45) rest frame and between 0.5 
and 5.0 GeV/c in the laboratory frame, and a polar angle 
in the range between 0.5 and 2.6 radians in the labo- 
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ratory frame. The total error on A z, computed on an 
event-by-event basis, is required to be less than 175 /im. 
The vertex fit constrains the lepton tracks to originate 
from the same point in the transverse plane, thereby ne¬ 
glecting the non-zero flight length for B mesons. As a 
consequence, the Az resolution function is A z depen¬ 
dent, becoming worse at higher \Az\. Neglecting this 
Az dependence introduces a small bias, discussed below. 
Vertex studies with leptonic J/ip decays show that the At 
resolution function for signal events can be appropriately 
modeled with the sum of three Gaussian distributions in 
both data and Monte Carlo simulation. 

The separation between direct leptons and background 
from cascade decays is achieved with a neural network 
that combines five discriminating variables for each event 
and provides two outputs, one for each lepton, chosen 
to vary between 0 for cascade leptons and 1 for direct 
leptons. The discriminating variables are the momenta 
of the two leptons, their opening angle, the total visible 
energy, and the missing momentum of the event. All 
variables are computed in the T(45) rest frame. The first 
two variables are very powerful in discriminating between 
direct and cascade leptons. The third efficiently removes 
direct-cascade lepton pairs from the same B decay and 
further reduces contributions from photon conversions. 
An optimization study based on the minimization of the 
total error on A rrid leads to the requirement that both 
neural network outputs be greater than 0.8. 

The numbers of selected on-resonance and off- 
resonance events are 99010 and 428, respectively. The 
combined effect of all requirements gives a direct dilepton 
purity and efficiency of about 83% and 9%, respectively, 
based on Monte Carlo simulation. Semileptonic B decays 
in the Monte Carlo simulation have been modeled sep¬ 
arately for each charm meson involved. A parametriza- 
tion of HQET form factors || is used for B —> D*lu, 
the Goity-Roberts model [H is used for B —> D^niu, 
while the ISGW2 model is used for B —> Div and 
B —> D**iv. The measured branching fractions for de¬ 
cays to D** and 7 r states are fixed to their world 

averages 0 and unmeasured processes have rates that 
are inferred on the basis of isospin arguments. Events 
from BB decays are grouped in three topologies, each of 
which is assigned its own PDF with different At depen¬ 
dence and tagging properties. 

Direct dilepton events are described by the convolution 
of an oscillatory term for neutral B decays, or an expo¬ 
nential function for charged B decays, with the resolution 
function 77: 


. |A t|/r„(c) 

5 »(c) = -- (1 ± Dj c) ^ (c) ) ® 77 


4t 

' r, 


(c) 


for neutral (n) and charged (c) events, where r„( c ) is the 
B meson lifetime, = cosAwdAf and £ c = 1, and 
£3 0.95 are correction factors that account for the 


(small) fraction of wrongly tagged direct dilepton events. 
These events are due to hadrons from the B vertex that 
are misidentified as leptons or leptons from the decay 
of resonances ( e.g ., events where only one lepton comes 
from a J/i/j) produced at the B vertex. Both of these 
sources give almost random tagging and, in the absence 
of such events, would be exactly 1. A small fraction 

of events where a lepton originates from the b —► —> i~ 
decay chain have the same charge correlation as signal 
events and are also included in the signal topology. Ne¬ 
glecting the r lepton lifetime introduces a negligible bias 
on the Arrid measurement. 

Opposite B cascade (OBC) events, 9% of the selected 
sample, contain one lepton from a b £ decay and the 
other from the b —► c —> £ decay chain of the companion 
B meson. These events are the main source of wrong 
tags. Their PDFs are modeled by the convolution of A t- 
dependent terms of a form similar to the signal with a 
resolution model that takes into account the effect of the 
charmed meson lifetimes by convoluting three Gaussians 
with a single-sided exponential decay distribution. Since 
both short-lived D° and D s , and long-lived D + mesons 
are involved in cascade decays, the global OBC PDFs are 
written as 


r n(c) _ 
'"'OBC — 


e -|At|/r n(c) 
1 A.(c) 


Y J fi {C \ l±D OVctn { c))®K OBC 


where the index i runs over the short- and long-lived 
charm meson components. This parameterization of 
OBC events significantly reduces the related systematic 
uncertainty. The two resolution functions 77 q BC allow for 
different effective charm lifetimes and parameters of the 
three Gaussians, since the resolution function depends 
on the B and D flight lengths. Due to the different 
decay processes involved, the relative fractions of 
short- and long-lived charm mesons are also different in 
neutral and charged B events. If particle identification 
were perfect and cascade leptons originated only from the 
b —> c —> £ + process, then flavor tagging would always 
be wrong and the factors D would be exactly — 1. 
Hadron misidentification (PID) and resonance decays, as 
well as leptons originating from the b —> cc(—> £~)s chain, 
give a fraction of right tags (15%) even in the OBC topol¬ 
ogy. These two processes have been factorized by writing 
Dobc = -Dpnj • an d assuming no correlation be¬ 

tween the two terms. 


Same B cascade (SBC) events, 4% of the selected sam¬ 
ple, contain two leptons from a single B meson, obtained 
via the decay chain b —> c£~V, with c —> SBC 

events are insensitive to mixing and, in the case of perfect 
particle identification and in the absence of resonances, 
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would always give opposite-sign leptons. The PDFs are 


r n(c) _ e 
'-'SBC — _ 


l At l/ T SBC ( C N 

(1 ± D sbc ), 


4r c 


»(c) 


SBC 


where are effective lifetimes and Dggg are correc¬ 
tions for wrong tags in the SBC topology. The resolution 
1Z is taken to be the same as for signal events, with no 
significant bias on the final result. 

A small residual background remains (0.3% of the total 
sample) where both leptons are from an unrecognized 
J/^ decay. These are described by a term T = 5{At)®TZ, 
whose normalization is obtained from simulation. Events 
where one lepton originates from a cascade decay and the 
other from a B decay to r or to a resonance, and events 
where both leptons come from cascade decays, (0.3% of 
the total sample) are assigned the OBC event topology 
with no significant bias on A rrid- 

The fraction / cont = 3.4% and At dependence of 
the continuum background are determined from off- 
resonance data. The At dependence is parameter¬ 
ized for opposite- and same-sign leptons as Q± = 
%on t e“ Tcont!At| /±, with /+ + /_= 1. 

The full likelihood function is the product of likeli¬ 
hoods for opposite- and same-sign events, which can be 
schematically written as 


£ = (1 /cont)(l fj/ip)[ 

(1 ~ /c)(/5gS" + /OBC^OBC + /SBC^SBC) + 
+W^S C + /obC^OBC + /sBC^SBc)] + 

T(1 ““ /cont) fjj Ik T /cont 

where the J/'i/j term and its relative abundance fj/^ 
are present for opposite-sign events only, and / s "^ = 

(1 — foBC ~ fsBG ) * The fraction f c of charged B events 
in the selected sample and the OBC fraction /q B q in 
neutral B events are extracted from the fit. The OBC 
fraction /q BC in charged B events is scaled with /q BC ac¬ 
cording to the value of the ratio /obc//obc determined 
with the Monte Carlo simulation. The SBC fractions 
are computed for simulated events and fixed in the fit. 
The various parameters for the OBC resolution functions 
are taken from a fit to Monte Carlo events. The factor 
Dp lB> (D°, D s ) is fitted and all the other corrections for 
wrong tags scale with D b1d (D °, D s ) according to ratios 
determined with simulated events. 

To summarize, the values for A md, fc , /obo 
D(1 ID (D 0 ,D S ), fa , and the widths and relative frac¬ 
tions of the Gaussian components for the signal resolu¬ 
tion are determined in the likelihood fit. The B meson 
lifetimes are fixed to the values quoted in [Jrj. 

The result of a binned maximum likelihood fit to the 
data sample with the requirement \At\ < 12 ps yields 
A m d = 0.488 ± 0.012 ps" 1 and f c = 0.554 ± 0.014. Fig¬ 
ure W and [lb show the At distributions for opposite- 



At (ps) 

FIG. 1: Distributions of decay time difference for (a) opposite- 
sign and (b) same-sign dilepton events; (c) asymmetry be¬ 
tween opposite- and same-sign dilepton events. Points are 
data and the lines correspond to the fit result. 


and same-sign dilepton events, respectively, along with 
the result of the fit. Figure [([• shows the resulting asym¬ 
metry as a function of At. The widths of the three Gaus- 
sians for the signal resolution function are 0.55 ± 0.09, 
1.06 ± 0.23 and 4.8 ± 0.7 ps, and the corresponding frac¬ 
tions of events are 76%, 22% and 2%. The probability to 
obtain a worse fit is 65%, evaluated with an ensemble of 
data-sized experiments that are generated with a param¬ 
eterized simulation based on the observed total PDF. The 
global fit is also performed on a sample from full Monte 
Carlo simulation, where the fitted results for parameters 
are consistent with generated values. 

The fit result is found to be stable and consistent under 
a variety of choices for free parameters, where fixed values 
obtained from Monte Carlo simulation are substituted. 

A summary of the systematic uncertainties is given in 
Table |, where the total is estimated to be 0.0087ps _1 . 
The most important contributions are due to the B me¬ 
son lifetimes, the At resolution function, and the mod¬ 
eling of OBC events. Varying the B meson lifetimes 
within their known errors [?]] contributes an uncertainty 
of 0.0064 ps -1 on A md- 

The systematic error due to the uncertain knowledge 
of the resolution function for OBC events is estimated by 
varying the parameters within their errors from the fit to 
simulated events, including the effect of correlations. A 
possible scale uncertainty between data and simulation is 
estimated by allowing a conservative increase of 20% in 
the OBC resolution width. The overall uncertainty due 




















7 


TABLE I: Summary of systematic uncertainties 


Source 

a(Am d ) ps 1 

B lifetimes 

0.0064 

OBC resolution/lifetimes 

0.0026 

At dependence of resolution 

0.0043 

z scale and SVT alignment 

0.0020 

OBC fractions/wrong tags 

0.0020 

Hadron misidentification 

0.0010 

J/ip fraction 

0.0003 

Continuum parameterization 

0.0009 

Binned fit bias 

0.0006 

Beam energy uncertainty 

0.0005 

Total 

0.0087 


giving a combined systematic error of 0.0009 ps” 1 . The 
dependence of the fit result on the number of bins has 
been estimated with a parameterized Monte Carlo simu¬ 
lation. A shift of —0.0006 ps -1 in Ais observed and a 
corresponding correction applied with a systematic error 
of 0.0006 ps -1 . The uncertainty (0.1%) on the absolute 
scale of the beam energies gives an error of 0.0005 ps” 1 
on A m.d- 

In conclusion, the neutral B meson oscillation fre¬ 
quency has been measured with an inclusive dilepton 
sample to be 

A m.d = 0.493 ± 0.012(staf) ± 0.009(sysf) ps” 1 . 


to the OBC resolution function is 0.0026 ps” 1 . The im¬ 
pact of our treatment for the At dependence of both the 
signal and OBC resolution functions and the boost ap¬ 
proximation has been studied with large parameterized 
Monte Carlo samples, which are based on the observed 
dependence in full simulation. Neglecting the At depen¬ 
dence results in a bias for A rrid of —0.0045 ps” 1 . The fit 
result has been corrected to account for this bias and a 
corresponding systematic error of 0.0043ps” 1 is assigned. 
Knowledge of the absolute 2 scale of the detector and the 
residual uncertainties in the SVT local alignment give a 
combined error of 0.0020 ps” 1 . 

Systematic effects due to the limited knowledge of the 
parameters of the OBC PDF, which are taken from sim¬ 
ulated events, are greatly reduced by fitting the frac¬ 
tions of OBC and the short-lived charm component in 
neutral B events. The remaining systematic uncertainty 
(0.0020 ps” 1 ) is estimated by varying the otherwise fixed 
charm-related parameters (the amount of D s , /^ 0 D and 
the various fractions of cascades) by 10%, both coher¬ 
ently and in independent directions. This is a conserva¬ 
tive range, given our present knowledge of the physics 
processes involved. 

The ratios between the various wrong-tag factors due 
to hadron misidentification (PID) are conservatively var¬ 
ied by 30% in the fit. The maximum effect is obtained 
when the signal and cascade PID wrong-tag corrections 
are varied in opposite directions. In this case, the total 
systematic error is 0.0010 ps” 1 . 

The uncertainty on the fraction of J/ip is 30%, which 
contributes an error on A m d of 0.0003 ps” 1 . The effec¬ 
tive lifetime, the fraction of same-sign events, and the 
fraction of continuum events are varied independently, 


This result is the single most precise measurement to date 
and is consistent with a recent BABAR measurement with 
a fully reconstructed B° sample || . 
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